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Abstract. We investigate a model of epidemic spreading with partial immunization 
which is controlled by two probabilities, namely, for first infections, po, and reinfections, 
p. When the two probabilities are equal, the model reduces to directed percolation, 
while for perfect immunization one obtains the general epidemic process belonging to 
the universality class of dynamical percolation. We focus on the critical behavior in 
the vicinity of the directed percolation point, especially in high dimensions d > 2. It 
is argued that the clusters of immune sites are compact for d < 4. This observation 
implies that a recently introduced scaling argument, suggesting a stretched exponential 
decay of the survival probability for p = pc, po ^ Pc in one spatial dimension, where Pc 
denotes the critical threshold for directed percolation, should apply in any dimension 
d < 3 and maybe for d = 4 as well. Moreover, we show that the phase transition line, 
connecting the critical points of directed percolation and of dynamical percolation, 
terminates in the critical point of directed percolation with vanishing slope for d < 4 
and with finite slope for d > 4. Furthermore, an exponent is identified for the temporal 
correlation length for the case of p — Pc and po — pc — e, e I, which is different 
from the exponent of directed percolation. We also improve numerical estimates 
of several critical parameters and exponents, especially for dynamical percolation in 
d = 4,5. 
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1. Introduction 

The study of stochastic models that describe the spreading of a nonconserved agent is 
currently a very active field of research in nonequilibrium statistical physics Pl El Ej- 
Fundamental interest stems from the fact that these models exhibit continuous 
nonequilibrium phase transitions from a fluctuating active phase into one or several 
absorbing states. Moreover, such models are motivated by a variety of possible 
applications [4, such as epidemic spreading jSl Ej, catalytic reactions [2,, or flowing 
sand HEl. 

Usually models for epidemic spreading are defined on a lattice whose sites can 
be either occupied (active, infected) or empty (inactive, healthy). The dynamic rules 
involve two competing processes, namely, spreading of activity to neighboring sites 
(infection) and spontaneous decay (recovery). Depending on the relative frequency 
of infection and recovery the process either has a finite probability to survive in an 
infinite system or it will eventually die out. Since the spreading agent is not allowed to 
be created spontaneously, once the empty lattice is reached, the process is trapped in 
a so-called absorbing state from which it cannot escape. Usually the transition between 
infinite spreading and recovery is a critical phenomenon characterized by large-scale 
fluctuations associated with certain universality classes. Therefore, the classification of 
all possible types of phase transitions into absorbing states is presently one of the major 
goals of nonequilibrium statistical physics. 

The most prominent universality class of phase transitions into an absorbing state is 
that of directed percolation (DP) |Tn|, as described by Reggeon field theory [TT j IT2 | IT8] . 
Models in this class describe short-range spreading of a non-conserved agent in a 
medium without temporal memory effects. Additionally, it is assumed that the rates 
for spreading and recovery are constant in space and time (no quenched randomness). 

If the medium has a memory such that each site can be infected only once, the 
transition does no longer belong to the DP class, instead one obtains the so-called 
"general epidemic process" (GEP) [SI IE!- In the language of epidemic processes this 
type of memory accounts for perfect immunization. In contrast to DP, where sites can be 
reinfected without restriction, a general epidemic process can only spread in those parts 
of the system which have not been infected before. Nevertheless, in spatial dimensions 
d > 2, infinite spreading from a single infected site in a non- immune environment is 
still possible, provided that the susceptibility to primary infections is sufficiently large. 
In this case activity propagates as a front, leaving a cluster of immune sites behind. 
The transition between survival and extinction of the spreading agent is then described 
by a different type of critical behavior, belonging to the universality class of dynamical 
percolation (DyP) ^Hl- In fact, the resulting cluster of immune sites has the same (or 
even identical) structure as an isotropic percolation cluster. For this reason DyP can be 
used as a dynamical procedure to generate isotropic percolation clusters. 

As an immediate generalization one may consider an epidemic process in which 
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Figure 1. Phase diagram of the epidemic process with finite immunization on a 
square lattice in 2+f dimensions, po denotes the primary infection probabihty and p 
the reinfection probabihty. The inset shows the vicinity of the DP point. 



the strength of immunization can be varied JHl- For example, the initial susceptibility 
to infections po may be locally set to a different value p when the first infection is 
encountered. In the following we shall denote this process as epidemic process with 
finite immunization which we abbreviate for convenience as EPFI. The phase diagram 
of the EPFI in two spatial dimensions was studied in jU^- As shown in Fig. ^ it 
comprises three different phases, including the GEP {p = 0) and DP {p = po) as special 
cases. 

The horizontal phase transition line in Fig.^can be explained as follows. Obviously, 
when starting from a fully occupied lattice, all sites become immediately immune so that 
the infection rate is everywhere equal to p. Trivially, the dynamics is then precisely that 
of a DP process controlled by p. This gives rise to a DP transition dXp = Pc, independent 
of pq, where pc denotes the critical value of an ordinary DP process. However, starting 
with a localized infected seed in a non-immune environment, the situation is more subtle, 
though the horizontal line still exists. For example, slightly above the critical line, i.e. 
for p = Pc + e with e -C 1, the interior of a surviving cluster is essentially dominated by 
an ordinary supercritical DP process. Since such a DP process is characterized on large 
scales by a finite density of active sites, there will be a finite chance for the process to 
survive in the limit t ^ oo |17j . 

The situation below the horizontal line {p < pc) is different. In at least two or 
more dimensions this part of the phase diagram displays two distinct phases, namely, an 
absorbing phase, where the process stops after some time, and a phase of annular growth. 
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where an expanding front of high activity survives with finite probabihty. Performing 
a field-theoretic renormahzation group study close to the upper critical dimension of 
DyP, dc^^ = 6, and computing the corresponding critical exponents, Janssen showed 
that the critical behavior along the line connecting the GEP and the DP point in Fig. ^ 
is that of DyP ^H] • In this sense immunization is a relevant perturbation, driving the 
system away from the DP point in Fig. ^ towards GEP. (We note that in the case 
of several competing infections with immunization, one observes a crossover back to 
DP [Hj. Another four-state generalization leads to in interesting tricritical phenomenon 
belonging to a different universality class [IH])- 

More recently epidemic spreading with immunization appeared in a different 
context. Studying systems with infinitely many absorbing states such as the non- 
diffusive pair contact process [201 1211 , Muiioz et al conjectured that these models are 
described by the same type of field theory as the EPFI |22l IIB] • Roughly speaking, 
the frozen absorbing configurations generated by the process provide a local memory 
of activity in the past which effectively acts in the same way as immunization (or 
weakening). Meanwhile this conjecture is widely accepted, although some questions 
concerning the upper critical dimension are still debated flM I25j . 

While the mentioned field-theoretic results explain the transition line between 
the GEP and the DP point close to d^^^ = 6, the critical properties along the 
horizontal line and in the vicinity of the DP point are less well understood. Simulating 
the Langevin equation of the EPFI at critical reinfection rate (corresponding to the 
horizontal line in Fig. in 1+1 dimensions starting with a localized seed, Lopez and 
Munoz initially expected continuously varying exponents [22], but refined simulations 
and approximations suggest that there is no power-law scaling. Instead, the activity 
was found to vary as a stretched exponential in 1+1 dimensions [T7| 12 7j. 

In this paper we study the epidemic process with finite immunization, focussing 
on the influence of immunization in the vicinity of the DP point, especially in higher 
dimensions d > 2. In Sec. El we describe Monte Carlo simulations that are applied in 
the following sections to obtain numerical results. In Sec. El we summarize the most 
important phenomenological features of the process. In particular, we argue that the 
clusters of immune sites at the DP point are compact for d < 3 and asymptotically 
compact for d = 4 (see below). Consequently, a recently introduced scaling argument 
that suggested a stretched exponential decay on the horizontal line for po <^ pc in 
d = 1 [27\ should apply in any dimension d < 3 and maybe for d = 4 as well, since it is 
based essentially on the compactness of the cluster of immune sites. This means that we 
expect a stretched exponential decay instead of power-law behavior not only in one but 
also in two, three and perhaps four spatial dimensions. Moreover, for d < 4, primary 
infections and reinfections scale differently whereas they show the same scaling behavior 
for d > 4. As a result, we find that the phase transition line connecting the GEP and the 
DP point should terminate in the DP point with a finite slope for c? > 4 (in contrast to 
a vanishing slope for d < 4, as can be seen in Fig. Q). In Sec. IHwe discuss the Langevin 
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equation for EPFI in order to describe how the properties of the process change between 
the upper critical dimension of DP, d^^ = 4, and DyP, dc^^ = 6. Numerically determined 
phase diagrams for d > 2 are presented in Sec. El which support the predictions of Sec. 01 
Though we expect the behavior along the horizontal phase transition line to be non- 
universal (at least for o? < 4), it is nevertheless possible to identify an exponent fi for the 
temporal correlation length in the case of p = Pc and po = Pc — e for e <^ 1. As discussed 
in Sec. El this exponent differs from the critical exponent for the temporal correlation 
length of DP. The main results are summarized in Sec. [3 A heuristic argument for the 



value of /i that is in accordance with the numerical data is presented in Appendix A 



Numerical estimates for critical parameters and critical exponents for DP and DyP are 
given in Appendix B and Appendix C[ respectively. 



2. Monte Carlo simulations 



We perform Monte Carlo (MC) simulations of the epidemic process with finite 
immunization on a simple d-dimensional cubic lattice using the model of Ref. ^7] 
generalized to arbitrary dimensions. The process is defined such that it reduces to 
directed bond percolation at the DP point and to dynamical bond percolation at the 
GEP point. In particular, an active lattice site i at time t may activate each of its d 
nearest neighbors, denoted as j, at the next discrete time step t + 1. If neighbor j has 
never been active before, activity is transmitted from site i to site j with the primary 
infection probability po. Contrarily, if site j has been active at least once in the past the 
infection probability is given by p, called reinfection probability. Each site stays active 
only for one time step. 

The process is initialized with an active seed at the origin at time t = in a non- 
immune environment. Each run is stopped either when the process dies out or when it 
reaches a preset maximum time. We average over many runs with different realizations 
of randomness. The lattice is always chosen large enough so that the process never 
reaches its boundary. Hence, finite size effects are eliminated. As usual in such 'seed 
simulations' we measure the survival probability Ps{t) that the process survives at least 
up to time t, the number of active sites averaged over all runs N{t), and the mean 
square spreading from the origin R'^{t) averaged over surviving runs. At criticality these 
quantities obey power-laws j2Hl 1211 according to 

P,(t)~t-^, N{t)^t' , R\t)^t^/' . (1) 

This behavior is observed for critical DP as well as for critical DyP. Of course, the 
exponents in Eq. are generally different in both cases. In addition, we also measure 
the number of primary infections (see below). 

From a technical point of view, our simulations for d > 3 are based on the routine 
presented in |^. Lattice sites are labeled by 64-bit-long integers. We do not initialize 
storage for a whole L'^ lattice since for high dimensions it is only possible to simulate 
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Figure 2. Typical cluster of immune sites generated by the EPFI in 2+1 dimensions 
at the DP point after 2048 updates. The colors represent linearly the time at which 
the sites were visited for the first time. The black line marks the final boundary of the 
cluster. 

small lateral lattice sizes L with this method. Instead we use lists to store the individual 
positions of active and immune sites. To perform an update from t to t + 1 we go through 
all active sites at time t and activate their neighbors at time t + 1 with probability po 
or p, respectively. The activated sites are stored in a list while the formerly active sites 
recover. In the case of a first infection, the site is added to the list of immune sites. In 
order to efficiently check whether a site is immune or not, and whether it has already 
been activated during the actual update step, we use a hashing algorithm as described 
in ISO]- Contrarily, the simulations for d = 1,2 are carried out in the usual way where 
one initializes storage for a whole lattice. However, in this case we apply bit-coding, 
i.e., a single 64-bit-long integer stores the states of 64 lattice sites (active/inactive or 
immune / non- immune) . 

3. Phenomenological properties 

Let us first consider the phenomenological properties of the process near the DP point, 
where the effect of immunization is small. Since the DP point itself is a critical DP 
process, the question arises how the upper critical dimension of DP influences the 
generated cluster of immune sites. As a conjecture we propose that the generated 
cluster of immune sites at the DP point is compact in d < 4 dimensions (in the sense 
that its fractal dimension is d), while in higher dimensions > 4 it is not. At the DP 
point p = Po = Pc, the influence of immunization vanishes so that the cluster of immune 
sites is merely given by the sites visited in the past by a critical DP process. In other 
words, when projecting the past activity of a critical DP cluster onto space (looking 
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Figure 3. Double logarithmic plot of the number of primary infections versus time at 
the DP point for various dimensions. For comparison the power-law behavior predicted 
by Eq. (PJ is shown by the dashed lines for d = 2,3. For d — i the mean-field prediction 
of Eq. JSJ is np{t) = const. 



through a critical DP cluster along the temporal axis) we expect its appearance to be 
compact in d < 4. Obviously, in = 1 the cluster of visited sites is compact by definition 
so that the conjecture is correct. In dimensions d = 2,3 the statement is non-trivial. 
In fact, plotting the projection of a typical cluster in d = 2 dimensions one obtains a 
compact object, as shown in Fig. |2l 

The above conjecture can be supported numerically as follows. Assuming 
compactness, the number of immune sites in a surviving run should grow linearly with 
the volume ^f, where ~ t^^^ is the spatial correlation length and z = i^\\/i^± is the 
dynamical exponent of DP (see, e.g., [2]). Averaging over all runs, the volume has to be 



multiplied with the survival probability Ps{t) 
number of immune sites I{t) increases as 



t , where 6 = (3/v\\, hence the average 



i{t) ~ r/"-' . (2) 

This implies that the number of primary infections scales as the derivative of I{t), 
i.e., 

^d/z~S-l 



np{t) 



(3) 



As shown in Fig. IHl by numerical simulations, in d = 2,3 this quantity scales indeed 
according to Eq. Q, supporting that the generated clusters of immune sites are compact. 



For d = A the data still exhibits a slight curvature which is presumably due to 
logarithmic corrections [211 IS21 that affect the mean-field behavior {up = const) at the 
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upper critical dimension of DP, d^^ = 4. As usual, also for d > d°^ we observe mean-field 
behavior and thus we expect the resulting cluster of immune sites to be non-compact. 

It is also possible to support the conjecture by a scahng argument. Obviously I(t) 
cannot be larger than the integrated past activity /q dt'N{t') ~ ^HI, leading to the 
inequality 

d/z-5<e + l . (4) 

as a necessary condition for compactness. In ci < 4 dimensions the initial-slip exponent 

is given by the hyperscaling relation ^28, 3Jlj 

e = d/z- 26 (5) 

so that the inequality reduces to 5 < 1, which is indeed satisfied in c? < 4 (see Tab. IB Hi . 
At the upper critical dimension of DP, d^^ = 4, taking the mean-field exponents ^ = 0, 
z = 2, and 6 = 1, the inequality (jlj) is sharply satisfied while it is violated above four 
dimensions. Hence, we can conclude that in c? > 4 dimensions the generated cluster is 
no longer compact. 

To deal with the question whether clusters of immune sites in systems at the upper 
critical dimension of DP, d^^ = 4, are still compact, it is worthwhile to compare the 
situation with a simple random walk whose upper critical dimension is d^^ = 2. The 
random walk is recurrent in d = 1,2 while it is transient for d > 3 j^l]. This means that 
the probability F{t) to visit a given site at least once until time t tends to 1 in the limit 
t ^ oo for d = 1,2 while it tends to a constant c < 1 for d > 3. Therefore, the region 
visited by the random walker is compact in d < 2 whereas it is not in > 3. However, 
ioT d = 2 the probabihty 1 — F(t) that the site has not yet been visited decreases 
asymptotically as 1 — F{t) ~ (lnt)~^ in contrast to algebraic behavior in d = 1, where 

1 — F{t) ~ This implies that in = 2 compactness is reached only slowly and 
we refer to this as asymptotically compact. Therefore, in ci = 2 and finite time the 
cluster of visited sites is highly non-compact at its boundaries and differs significantly 
from the one shown in Fig. El Using this analogy, it is near at hand to speculate that 
the same happens in a critical DP process. Thus, we expect the cluster of visited sites 
to be compact at the upper critical dimension d^^ = 4 as well. However, in this case 
compactness should be reached considerably slower than for d < 4, i.e., we expect the 
clusters to be asymptotically compact. 

The observed compactness in low dimensions leads to an important consequence 
regarding the critical behavior along the horizontal line in Fig. ^ as it allows the results 
of ^22j obtained in d = 1 to be generalized to 2 and 3 dimensions. In |23 the expansion 
of an immune region in one spatial dimension was studied in the limit of a very small 
probability of first infections po ^ 0. Using a quasi-static approximation it was shown 
that the survival probability does not obey a power law, instead it decays as a stretched 
exponential 

Ps{t) oc exp (-Apo "t^-") , (6) 
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giving rise to a finite survival time T ~ Pq^^^ "\ Here a = I'w/ii'x+Ps) with f3s being the 
order parameter exponent next to a planar absorbing surface jHSj- The main assumption 
made in this approximation is the compactness of the immune domain, wherefore the 
analysis had been restricted to d = 1. Because of the observed compactness, the same 
arguments can be applied in 2 and 3 dimensions so that the formula Eq. (jH)) should be 
valid in these cases as well. Inserting the known values for /3<j ^| one obtains a ~ 0.947 
in Id and a ~ 0.72 in 2d (so far there are no estimates of in 3(i). Whether Eq. © 
is still valid at the upper critical dimension of DP, d^^ = 4, where a would be 1/2, is 
not yet clear since in this case we expect the immune region to be only asymptotically 
compact. The numerical verification of the stretched exponential behavior for d = 1 and 
Po ^ Pc in [2Z| required a special enrichment method since for small po most simulation 
runs terminate after a very short time. In this paper we perform ordinary Monte Carlo 
simulations and therefore we do not address the numerical analysis of the horizontal line 
for Po <^ Pc in d > 1. 

Another important consequence concerns the ratio of first infections and 
reinfections. In < 4 dimensions the average number of first infections np{t) (controlled 
by the parameter po) decreases according to Eq. while the average number of 
reinfections N{t) (controlled by the parameter p) increases as with > 0. Thus, 
after sufficiently long time the parameter p controlling reinfections will have a much 
larger influence than the parameter p^. As a consequence the curved phase transition 
line, which may be thought of as describing a situation where p and po balance each 
other, terminates in the DP point horizontally, i.e. with zero slope. In > 4 dimensions, 
primary infections and reinfections show the same scaling behavior (for d = A probably 
affected by logarithmic corrections), and therefore we expect the transition line to 
terminate with a non- vanishing slope. We will come back to this question in Sec. IHl 



4. Scaling properties 



The previous phenomenological arguments suggested that the behavior of the model 
in d < 4 dimensions differs significantly from the behavior above 4 dimensions. This 
difference becomes also obvious when studying the corresponding Langevin equations. 

According to Refs. ^Hl HZj the Langevin equation for the epidemic process with 
finite immunization reads 

^p(x, t) = ap(x, t) - 6p2(x, t) + Z^Vp(x, t) + e(x, t) (7) 

+ Ap(x, t) exp(^-wj^ dTp{-x,T)^ , 

where ^ (x, t) represents a density-dependent Gaussian noise with the correlations 

(^(x, t)e(x', t')) = F p(x, t) 5^(x - x') 6{t - t'). (8) 

It consists of the usual Langevin equation for DP plus an exponential term describing 
the effect of immunization. Here the exponential function can be thought of as a switch: 
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Initially, the integral is zero and hence the coefficient of the linear contributions in 
p(x, t) is a + A, representing as bare parameters the reduced rate for ffist infections 
Po — Pc- However, when the integrated past activity exceeds 1/w, the exponential 
function decreases rapidly so that the additional term is essentially switched off. Roughly 
speaking, in a continuous description the parameter w is needed in order to specify a 
threshold telling us how much activity has to be accumulated at a given site in order 
to declare it as immune. Once the exponential term is switched off, the linear term is 
controlled by the coefficient a which represents the reduced reinfection rate p — pc- In 
most lattice models, sites become immune after a single infection, hence the parameter 
w is of the order 1 while X ^ po — p controls the strength of immunization. 
Rescaling the DP Langevin equation by 

x^6x t^bH p^b-^'p (9) 

with a scaling parameter b and the exponents 2; = z/y / z/_l and x = l^/^i. immediately 
recognizes that simple scaling invariance at the upper critical dimension d^^ = 4 can 
only be established if 2; = x = 2 and a = 0, the latter representing the critical point 
at the mean field level. Regarding the exponential term, scaling invariance requires 
the argument of the exponential function and the exponential function itself to be 
dimensionless, hence the coefficient A and w have to be rescaled as 

A ^ b-y^X (10) 
w — > b^y^^w , 

where y\ = 2 and = 0. More generally, it can be shown by a field-theoretic 
calculation j22I that in ci < 4 dimensions the two exponents are given by 

1 z/|| - /3 
z/A = — , yu, = • (11) 

As Uw and y\ are positive, the exponential term for immunization is relevant m. d < A 
dimensions. Expanding the exponential function as a Taylor series, the resulting terms 
would be equally relevant in c? = 4 and increasingly relevant m d < A dimensions. 
Therefore, in < 4 dimensions a Taylor expansion of the exponential function in Eq. (jZj) 
is meaningless in the renormalization group sense, instead it has to be kept in as a whole. 
This circumstance is probably responsible for the observed non-universality along the 
horizontal line in Fig. ^ 

In (i > 4 dimensions, however, the situation is different. Here, power counting at 
the upper critical dimension of DyP, d^^^ = 6, yields y^ < 0, meaning that the relevancy 
of the terms in the Taylor expansion decreases. In this case it is legitimate to carry out 
the Taylor expansion, keeping only the most relevant contribution. As the zeroth order 
can always be absorbed in a redefinition of a, the most relevant contribution is the ffist- 
order term —Xwp{x, t) /q drpix, r). In addition we may drop the quadratic contribution 
—bp'^{'x,t), which is irrelevant at dc^^ = 6, leading to the Langevin equation 

^p(x, t) = {a + A)p(x, t) + DVp(x, t) + e(x, t) (12) 
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— Xw p{x,t) / dTp{-K,T). 
Jo 

corresponding to the field theory studied in [inillHl- In contrast to Eq. the influence 
of immunization is effectively described by a single parameter, namely, by the product 
Xw. Assuming to be negative for any d > 4 this scenario is expected to hold even 
in presence of fluctuation effects. Moreover, in contrast to Eq. ((7j) the parameter now 
appears in the prefactor of the linear term, shifting the parameter a. This suggests, 
in accordance with the preceding section, that at least in 4 < (i < 6 dimensions the 
reduced flrst infection and reinfection probability (here corresponding to A and a) scale 
identically and that the phase transition line, along which the influence of a and A is 
balanced, terminates with a nonzero slope in the DP point. 

To summarize we arrive at the following picture: 

• In d < 4 dimensions the EPFI is described by Eq. (j?)), in which an expansion of 
the exponential function is not allowed. The exponential function is conjectured to 
produce non-universal features such as a stretched exponential decay of the survival 
probability, as observed in d = 1 dimensions [27] . 

• In 4 < d < 6 the process is described by Eq. (fT^ . Fluctuation effects are still 
present and the corresponding fleld theory is well-deflned and renormalizable jl8j . 

• In d > 6 the mean field approximation becomes valid. The system is driven to a 
trivial Gaussian fixed point, which is the same for GEP and DP. 

5. Phase Diagrams 

Fig. m presents phase diagrams of the epidemic process with finite immunization for 
spatial dimensions d = 3,4,5,6 which were obtained obtained from MC simulations 
(see Sec. The horizontal lines obey p = Pc with the values of Pc given in Tab. IBll 
The lines connecting the DP and the GEP point were determined as follows. First we 
performed spreading simulations at the critical point of the GEP, i.e., for p = 0, using 
the critical values of bond percolation given in jSOj. Thereby, we estimated various 
critical exponents for DyP which are presented in Tab. 1(^11 In order to locate the phase 
transition line, we then used the fact that the critical behavior along this line is that of 
dynamical percolation. We determined various critical points along the lines by keeping 
p fixed and varying po until the quantities Ps(t), N(t) and R'^it), see Eq. (jT)), displayed 
the expected slope corresponding to DyP in a log-log plot. 

In accordance with the predictions of Sec. El the data indeed suggests that the 
curved phase transition line terminates at the DP point with vanishing slope for d = 3 
and with finite slope for d > 4. For d = 4 this behavior is not so clear cut and we believe 
that this is caused by logarithmic corrections at the DP point. 

Assuming that the curved phase transition line behaves in the vicinity of the DP 
point as po,c(p) — Pc ~ {Pc — pV , our data leads to 7 = 0.47(7) for ci = 2 and to 
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Figure 4. Phase diagrams for the EPFI for d — 3,4,5,6, A)-D), determined from 
MC simulations. The insets show the vicinity of the DP point. Circles mark the 
numerically determined critical points. 



7 = 0.74(8) for d = 3. Since the error bars are rather large, our data does only allow a 
rough estimate of 7 and we did not further attempt to relate its value to other critical 
exponents. 



6. Temporal correlation length in the vicinity of the DP point 

Off criticality but in the vicinity of the DP point, i.e., p = Pc^e, Pq = Pc^e, one initially 
observes the critical behavior of DP until the process eventually crosses over to a different 
type of behavior. This crossover takes place at a certain typical time scale ^\\. If e 7^ 
(moving vertically away from the DP point) one expects that asymptotically ~ e~'^ii 
with the critical DP exponent for the temporal correlation length. The reason is the 
following. Near the DP point, for d < 4, the process is dominated by reinfections that 
perform an ordinary DP process with parameter p on a compact region (see Sec. El). It is 
only on the edge of the visited region where the first infection probability po plays a role. 
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Figure 5. Data collapses for p = pc and po — Pc ~ £ for d — 2 (top) and c? = 3 
(bottom). N{t) is scaled with with the values of 6 from Tab. lBll The time is scaled 
with such that the best data collapses were obtained. The insets show the original 
data. 



Hence, the time scale associated with e is expected to be smaller than that associated 
with a change e of po- For d > 4, primary infections and reinfection show equal scaling 
behavior and therefore one expects both times to scale in the same way. 

We are left with the question how the temporal correlation length behaves for p = Pc 
and Pq = Pc ^ € (moving horizontally away from the DP point) in ci < 4. Although we 
expect non-universal behavior on the horizontal line it may nevertheless be possible to 
identify an exponent /i for the temporal correlation length 

eil ~ e-'' (13) 

in the vicinity of the DP point, e ^ 1. In fact, as shown in Fig. plotting N{t)t~^ 
versus t/^\\ we can produce a data collapse by tuning /x. With similar simulations in 
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d = 4,5 (not shown here) we get the estimates 



/i 



' 2.7(7) 


in d = 


2 


1.4(2) 


in d = 


3 


1.1(13) 


in d = 


4 


[ 1.0(1) 


in d = 


5 



(14) 



Obviously, in d < 4 the exponent fx differs from z/y (z/y = 1.295 for d = 2 [261 and 
z/|| = 1.105 [HZI for d = 3). Because of logarithmic corrections, it is likely that yU = 1 in 
d = 4 as well. With decreasing dimension the value of /i and therewith the simulation 
time increases rapidly. For this reason the estimate for d = 2 is less precise, while 
simulations in d = 1 turned out to be unreliable. However, for d = 3 the data collapse 
is fairly good. In App. [Appendix A we present a heuristic argument for the value of /i 
in accordance with the numerical data. 



7. Summary 

We studied a model for spreading with finite immunization which is controlled by 
probabilities for first infections (po) and reinfections (p). We focussed in particular on the 
critical behavior close to the directed percolation point, especially in high dimensions 
d > 2. We argued that we expect the domains of immune sites to be compact for 
d < 4, however, with an approach to compactness that is (logarithmically) slow at 
the upper critical dimension of DP, d°^ = 4. Therefore, we denoted the visited region 
as asymptotically compact in d = 4. The compactness of the immune region was 
supported by MC-simulations. We pointed out that this compactness implies that a 
recently introduced scaling argument, suggesting a stretched exponential decay of the 
survival probability for p = pc, po -C Pc in one spatial dimension, should apply in any 
dimension d < 3 and maybe in c? = 4 as well. Furthermore, we showed that for d < 4 
the number of first infections (averaged over all runs) decreases whereas the number of 
reinfections increases. Contrarily, for d > 4 both quantities scale equally. From that we 
derived the result that the phase transition line connecting the GEP and the DP point 
terminates in the DP point with a finite slope for d > 4 and with a vanishing slope 
for d < 4, which was supported by numerically determined phase diagrams. We also 
discussed the Langevin equation for the process, in order to study how the properties of 
the process change depending on the spatial dimension. Investigating the behavior for 
p = Pc and Po = Pc — e, e <^ 1 we were able to identify an exponent /x for the temporal 
correlation length which is different from z/y for DP. 
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Appendix A. Speculation about the value of /i 



In the following we present a heuristic argument for the value of /i in c? < 4 dimensions 
that is in accordance with the numerical results presented above. Consider first a 
subcritical DP process governed by the parameter p = e -C 1. Initially the system 

behaves as if it was critical until the hostile conditions eventually lead to extinction. This 
happens on a time scale ~ e~'^K During the active time the process produces active 
sites according to N(t) ~ t^, see Eq. (jT}. Therefore, a subcritical DP process activates 
on the whole M sites before it dies out, where M scales as M ~ Consider now 

the EPFI with p = Pc and Po = Pc ~ ^ which dies out on a time scale ^\\ ~ e^^. In this 
case (for d < A) one has basically a critical DP process in the compact interior region 
of the immune cluster and a subcritical process on its edge. The integrated activity on 
the edge scales as I{t) ~ f^^^'^, see Eq. (j2I). If one assumes that the critical reinfections 
do not introduce a time scale and that the process on the edge is subject to the same 
bound as subcritical DP concerning the maximal activity, one arrives at / ~ M leading 
to 

(A.1) 

d — oz 

With the critical exponents of DP given in Sec. IHl and in Tab. IBll Eq. ()A.1|) results in 
fi ^ 2.335 for d = 2, fi = 1.437 for d = 3 and fi = 1 for d = A which is compatible 
with the numerical results. For d = 1, with Tab. IBll and z/|| = 1.7338 [HSj one obtains 
fi ^ 4.814. However, this argument needs to be substantiated and therefore Eq. ()A.1|) 
has to be taken with care. 



Appendix B. DP in high dimensions 



In this Appendix we report some numerical results on directed bond percolation in high 
dimensions. The values are partly more precise than previously reported estimates. 
Similar results for directed site percolation in = 3, 4, and 5 dimensions were recently 
reported in 132] ■ 

Table IbT] shows the results for directed bond percolation ranging from 1 to 7 spatial 
dimensions. Our estimates for d = 3 differ slightly from that presented previously by 
Dickman gUI, namely, 6 = 0.7263(11), 6 = 0.110(1) and z = 1.919(4). Though the error 
bars in |1U] are smaller than in our case, our data is incompatible with the values of 
S and z in [10]. We do not conclude with respect to this issue but note that further 
numerical effort in this context is desirable. The value 6 = 0.107(5) is smaller than 
Dickmans result 6 = 0.110(1) but larger than the field-theoretic estimate 6 = 0.098 jH] 
based on a two-loop calculation. 

Finally we note that our estimates in 3d are compatible with the DP hyperscaling 
relation IHH] where, inserting our estimates for z and 6, we obtain 

e = d/z- 26 = 0.107(10) (B.l) 
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d 


Pc 


6 


6 


z 


Ref. 


1 


0.644700185(5) 


0.159464(6) 


0.313686(8) 


1.580745(10) 


EH] 


2 


0.287338(3)t 


0.4505(10) 


0.2295(10) 


1.766(2) 


EH] 


3 


0.1774970(5) 


0.737(5) 


0.107(5) 


1.897(5) 


present work 


4 


0.1288557(5) 


1* 


0* 


2* 


present work 


5 


0.1016796(5) 


1 





2 


present work 


6 


0.0841997(14) 


1 





2 


present work 


7 


0.07195(5) 


1 





2 


present work 



Table Bl. Numerical results for directed bond percolation on a simple d-dimensional 
cubic lattice. For d = 4, * denotes that the mean-field behavior is subjected to 
logarithmic corrections, z — v^^/v^. f Taken from |39) . 



d 


Pc 


6 


6 


z 


Ref. 


2 


0.5t 


0.092 


0.586 


1.1295 




3 


0.2488125(25) 


0.346(6) 


0.488(7) 


1.375(5) 


present work 


4 


0.1601310(10) 


0.595(8) 


0.30(1) 


1.605(9) 


present work 


5 


0.1181718(3) 


0.806(12) 


0.134(10) 


1.815(10) 


present work 


6 


0.0942019(6) 


1* 


0* 


2* 





Table CI. Numerical results for dynamical bond percolation on a simple d- 
dimensional cubic lattice. For d — 6, * denotes that the mean-field behavior is 
subjected to logarithmic corrections, z = i^^^/iy^. ^ Taken from f5 . Uncertainties 
are in the last digit. For d = 1 the transition is shifted to Pc — 1- 



In c? > 4, however, the hyperscahng relation is violated. 
Appendix C. Dynamical percolation in high dimensions 

Here we present our numerical estimates for the values of critical thresholds and 
exponents for dynamical bond percolation on a d-dimensional cubic lattice. Note that, 
as far as we know, no numerical values for critical exponents for DyP in c? = 4, 5 have 
been published before. The values presented for d = 3 are compatible with precise 
estimates pubhshed previously, i.e., = 0.2488126(5) 6 = 0.345(4), 6 = 0.494(6) 
and 1/z = 0.728(2) where the latter values are from 

The value of 6' = 0.488(7) is incompatible with the value of ^ = 0.536 reported 
in ^3] based on previous estimates from [IH]. However, in this case 6 was not directly 
measured but obtained via scaling relations, which presumably is the reason for the 
considerable difference between the two estimates. 

The estimates given in Tab. 1(^1 are compatible with the hyperscahng relation for 
DyP inSl iSl- For = 3, 4, 5 we obtain 

f 0.49(2) for = 3 , 
e = d/z~25-l = I 0.30(3) for = 4 , (C.l) 
[ 0.14(4) for = 5 . 
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